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Abstract: Monthly projections of maximum temperature, relative hu¬ 
midity, precipitation, and wind speed were made based on the model of 
HadCM3 and the climatic change scenarios of IPCC SRES A2a and B2a 
for the future scenario periods of 2010-2039 (referred to as 2020s), 
2040-2069(referred to as 2050s), and 2070-2099(referred to as 2080s). 
The period 1961-1990 was chosen as the baseline period. The observed 
and projected weather data were downscaled using delta change methods 
and historical relationships between weather data, area burned, and the 
seasonal severity rating (SSR) code of the Canadian Fire Weather Index 
System were examined. The variations of area burned as influenced by 
climate change were assessed quantitative and qualitative for the study 
region, assuming that the fire regimes had the similar responses to the 
warming climate as during the 20th century. Our results indicated that a 
linear regression relationship existing between the historical area burned 
and the mean SSR values with regression coefficient in the significant 
range of 0.16 to 0.61. It was evident that the increased SSR values could 
result in more area burned; the area burned in the study region would 
have an increasing pattern during the 21st century under scenarios A2a 
and B2a scenarios and the area burned would be doubled. Also, the future 
area burned would have a strong seasonal pattern that more fires would 
occur in summer and autumn fire season, especially in summer. The area 
burned in summer fire season would increase by 1.5 times compared to 
that in the baseline period in 2080s under A2a scenarios. 
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Introduction 

Fire is one of the dominant disturbances in China’s boreal forest. 
Although China’s forest cover has grown from approximately 
12% to more than 20% through sustainable forest management 
over the past 30 years (He et al. 2007; He et al. 2008; Zhu. 2009), 
the forest area affected by fires increases continually at an 
alarmingly rate of 0.6 million hectares per year. Area burned can 
be strongly influenced by a complex set of variables, such as 
climate/weather, ignition, fuel and fire control activity, especially 
climate/weather (Flannigan et al. 1991; Wheaton.2001). Recently, 
the climate has been warming as a result of an increasing con¬ 
centration of greenhouse gases and aerosols in the global atmos¬ 
phere composition (Douville, 2002). The global average surface 
temperature has increased by 0.74°C during the last 100 years 
(IPCC et al. 2007). China’s temperature rise has kept the same 
pace with global warming (Ding et al.2006). Such warming is 
likely to have a significant impact on area burned. It is suggested 
that the observed increase in area burned of China due to climate 
warming, for example, the average annual area burned of the 
total forest area increased from 0.5% in 1990-1999 to 0.8% in 
the period of 2000-2009. However, this warming has not been 
uniform across the whole country. Most of the warming in China 
over the last 50 years has been in the Northern Hemisphere upper 
latitudes, where the Great Xing’an Mountains boreal forest re¬ 
sides (Fang et al. 2000; Jia et al. 2009). The most productive 
forests area located in the Great Xing’an Mountains in north¬ 
eastern China, unfortunately, this region is always associated 
with high forest fire danger rating. Hence, study on the future 
area burned trends in the Great Xing’an Mountains is important 
in advancing the understanding of the climate change on forest 
fire regimes of boreal ecosystem. 

In recent years, there are a few methods available to estimate 
future area burned, and the use of the 21st century climate simu¬ 
lations of General Circulation Models (GCMs) and the historical 
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relationships between observed area burned and associated 
weather and fire weather index seems to be an appropriate 
methodology to apply (Flannigan et al. 2005; Tymstra et al. 2007; 
Carvalho et al. 2009; Krawchuk et al.2009). For instance, Flan¬ 
nigan et al. (2005) estimated that area burned in Canada may 
increase by 74%—118% by the end of 21st century using histori¬ 
cal relationships between weather/fire danger and area burned in 
tandem with two GCMs and these estimates do not explicitly 
take into account any changes in vegetation, ignitions, fire sea¬ 
son length or human activity that may influence area burned. 
Tymstra et al. (2007) estimated the area burned based on the 
historical frequency of area burned using Canadian Forest Fire 
Weather Index (CFFWI) system class. Carvalho et al. (2008) 
performed regressions for the area burned and fire occurrence 
with the meteorological and FWI variables. Balshi et al. (2008) 
developed relationships between air temperature and fuel mois¬ 
ture codes derived from the CFFWI system to estimate annual 
area burned using Multivariate Adaptive Regression Spline ap¬ 
proach across Alaska and Canada. 

The objective of this study was to predict the area burned in 
Great Xing’an Mountains boreal forests by forest fire for the 21st 
century. To our knowledge, a number of the studies on future 
area burned as influenced by climate change are underway, 
however, nearly all these studies are conducted in North America, 
Europe and no previous studies have addressed this issue in 
China. 

Materials and methods 

Study area 

The study area covers approximately 8.35 million hectares of 
Great Xing’an Mountains, northeastern China (50°10'-53°33' N, 
121°12'-127°00' E). The study region comprises three counties 
(Mohe, Tahe, Huma) and four districts (Jiagedaqi, Songling, 
Xinlin, Huzhong) (Fig.l). The region has a typical continental 
climate, with long, cold winters and short, cool summers. The 
mean annual temperature is -2.8°C and total annual precipitation 
is 746 mm. Elevation varies gradually from approximately 300 to 
700 m. The dominant species in the region area are Larix gmelini, 
Betula platyphylla , Pinus sylvestris var. mongolica, etc. (Xu, 
1998). It is observed that the species of higher dominance in the 
landscape area are more prone to wildfires. 

Data collection 

Forest fires data 

Forest fire records were obtained from the local fire management 
agencies of the Great Xing’an Mountains. The fire records from 
1966 to 2008 include the date of ignition, location, ignition cause, 
fire size. There are two fire seasons in this region: the spring fire 
season (15 May to 15 June) and the autumn fire season (15 Sep¬ 
tember to 15 November). Forest fires occur throughout the fire 
seasons; however, recently they are becoming more frequent and 
burning large areas in summer months. Many studies have fo- 
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cused on forest fires in summer, and have considered the time 
period from 16 June to 14 September as the summer fire season 
(Zhang 2008; Zhang 2009; Zhang et al. 2010). In this paper, all 
fire records during the three fire seasons were included in the 
analysis for consistency with the meteorological records. 



Fig. 1 Location and administrative division of Great Xing’an Moun¬ 
tains, northeastern China. 

Meteorological data 

The meteorological data for the study were obtained from China 
Meteorological Data Sharing Service System 
(http://cdc.cma.gov.cn/) for the same period as the fire data from 
1966 through 2008. Eighteen meteorological stations were se¬ 
lected for analysis (see Table 1). For each station, the data con¬ 
tained daily maximum temperatures, 24-h precipitation, mean 
wind speed, and mean relative humidity. Not all the stations have 
equally complete meteorological records (Table 1). To fill these 
gaps, the stations with intermittent meteorological data were 
estimated by linear regression from the nearest stations. 

Climatic scenarios 

The climate change simulations produced by the third version of 
the UK Meteorological Office’s Hadley Centre Coupled Model 
(HadCM3) were used to estimate the discharge of the study area. 
HadCM3 is the tool most commonly used to project future cli¬ 
mate and exhibits better performance in East Asia. The model 
requires no artificial flux adjustments to prevent excessive cli¬ 
mate drift and its dynamics and physics are solved on a 
3.75°x2.5° longitude-latitude model grid with 19 hybrid vertical 
levels (Carson 1999 and Lynch 2008). The climate change sce¬ 
narios were obtained from the Intergovernmental Panel on Cli¬ 
mate Change (IPCC) Data Distribution Centre 
(http://www.ipcc-data.org/). Within the IPCC Special Report on 
Emission Scenarios (SRES), the SRES A2a and B2a scenarios 
were used in this analysis. Both A2a and B2a describe a “region¬ 
alization”, meaning a heterogeneous world development. A2a 
describe a highly heterogeneous future world with regionally 
oriented economies. The B2a is also regionally oriented but with 
a general evolution towards environmental protection and social 
equity. Compared to A2a, B2a has a lower rate of population 
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growth, a smaller increase in GDP but more diverse technologi¬ 
cal changes and slower land-use changes (Nakicenovic et al. 
2000; Leckebusch et al. 2004). The variables available for 
HadCM3 were expressed as 30-year monthly mean values, cal¬ 
culated for the baseline period 1961-1990 (refers to here as 
1980s), and predicted fields for the 2010-2039 (refer to as 
2020s), 2040-2069 (refer to as 2050s), and 2070-2099 (refers to 
as 2080s) (Arnell 2003). The climatic variables available for the 
two scenarios include mean maximum temperature, mean pre¬ 
cipitation, wind speed and relative humidity. 


Table 1. The information of meteorological stations used in the sta¬ 
tistical analyses 


Station name 

Latitude 

Longitude 

Elevation(m) 

Start year 

Missing Year 

Mohe 

52.97 

122.52 

433.0 

1957 

1958 

Tahe 

52.35 

124.72 

361.9 

1960 

1962-1971 

Xinlin 

51.70 

124.33 

494.6 

1972 


Huma 

51.72 

126.65 

177.4 

1954 


Amuer 

52.83 

123.18 

550.0 

1974 


Huzhong 

52.05 

123.67 

520.8 

1974 


Eergunayouqi 

50.25 

120.18 

581.4 

1957 


Tulihe 

50.48 

121.68 

732.6 

1957 


Daxinganling 

50.40 

124.12 

371.7 

1966 


Heihe 

50.25 

127.45 

166.4 

1959 


Xiaoergou 

49.20 

123.72 

286.1 

1957 


Nenjiang 

49.17 

125.23 

242.2 

1951 


Sunwu 

49.43 

127.35 

234.5 

1954 


Boketu 

48.77 

121.92 

739.7 

1951 


Zhalantun 

48.00 

122.73 

306.5 

1952 


Beian 

48.28 

126.52 

269.7 

1958 


Keshan 

48.05 

125.88 

234.6 

1951 


Fuyu 

47.80 

124.48 

162.7 

1956 



Statistical analyses 
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are the mean monthly HadCM3 simulated future temperature(°C), 
precipitation(mm), relative humidity(%) and wind 


speed(m/s); T 


p,monthly 5 p,monthly 5 p,monthly 


and W 


p,monthly 


are the 


mean monthly HadCM3 simulated present temperature(°C), pre- 
cipitation(mm), relative humidity (%) and wind speed (m/s). 


Forest fire weather index 

The CFFWI system consists of six components: three fuel mois¬ 
ture codes (Fine Fuel Moisture code, Duff Moisture Code, 
Drought Code) and three fire behaviour indices (Initial Spread 
Index, Build Up Index, Fire Weather Index). The Fire Weather 
Index (FWI), a daily measure of potential fire danger, is derived 
from daily temperature, relative humidity, wind speed, and 24-h 
precipitation. The FWI is used to predict the intensity of a 
spreading fire, largely by estimating the state of forest fuels. The 
FWI is also used to estimate the fire hazard warning that is 
commonly posted in the Canadian northern communities. The 
Daily Severity Rating (DSR) is calculated from the FWI accord¬ 
ing to the following equation: 


Delta change methods 

Daily rather than monthly data were used because the weather 
and, consequently, fire behavior can change dramatically over 
time periods much shorter than a month. Downscaling methods 
were thus needed to bridge the gap between the large scale cli¬ 
mate scenarios and the finer local climate scale and to change 
monthly climate variability into daily local climate variability. 
Delta change approach was used in this paper, which is very 
simple and has been used in many climate change impact studies 
(Lenderink et al. 2007; Zhao FF et al. 2007; Akhtar et al. 2008). 
The monthly climatic variables (temperature, humidity, precipi¬ 
tation and wind speed), obtained for the period 1980s and pre¬ 
dicted mean monthly change fields for 2020s, 2050s, and 2080s 
for both A2a and B2a simulations, were downscaled by delta 
change approach for the regional model simulations. The four 
daily climatic variables for the future time series between the two 
simulations were computed according to the following equations: 

T —T +(T —T ) (1) 

f,daily o,daily \ f,monthly p,monthly / V / 


DSR = 0.0272 (FWI) 177 (5) 

The purpose of the DSR is to weigh the FWI to improve esti¬ 
mates of the effort required to control a fire under different con¬ 
ditions. By averaging DSRs over a period, one can obtain the 
Seasonal Severity Rating (SSR). The SSR, a final component of 
CFFWI system (Van Wagner C E, 1987), was used to measure 
fire seasonal variation within the climate of the region in this 
study. The SSR was calculated as the mean seasonal Daily Se¬ 
verity Rating of all meteorological stations in the region over the 
fire season. 

Results and discussion 

Area burned trends 

Historically fire has a strong impact on the vegetation in Great 
Xing’an Mountains. Current estimates indicated that an average 
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of 0.15 M ha was burned annually across the entire study region 
since 1966. Fires in 1966, 1972, 1977, 1987 and 2003 caused 
enormous damage to the area burned (more than 0.7 M ha). Al¬ 
though there is great year-to-year variability of average annual 
area burned, there is evidence that the area burned has decreased 
dramatically in recent decades as a result of strict enforcement of 
the fire suppression policy. The average area burned has de¬ 
creased more than 80% since 1990. Analysis of average areas 
burned by decades revealed a bimodal distribution with distinct 
low and extreme high area burned years (see Fig. 1). Statistical 
evidence suggested that there has been a decreasing trend in area 
burned starting in the early 1980-1989 period and continuing 
into the 1990-1999 period. The projected fire statistics showed 
an increase pattern from 2000 and onwards. 



Fig. 1 Annual mean area burned in the period of 1966-2008 for 
Great Xing’an Mountains. 

Based on the fire information provided, area burned in the 
study region is strongly seasonal, with 83.9% of area burned 
occurred in the spring fire season during 1966-2008. Also in the 
same period, a significant increase of area burned happened in 
the summer time. Since 2000, wildfires have burned approxi¬ 
mately 35 520 ha in the summer fire season around this region, 
which are accounted for 36.7% of that since 1966. The monthly 
time span of area burned ranged from the beginning of March to 
the end of October (see Fig.2). On average, May is the month 
with the highest area burned. 



month 


Fig. 2 Monthly mean area burned in the period of 1966-2008 for 
Great Xing’an Mountains. 
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Relationship between forest fires and fire weather indices 

To examine the relationship between area burned and the mean 
seasonal SSR value for the period 1961-1990, Linear regression 
analyses were conducted between area burned and the SSR val¬ 
ues with the latter as independent variable for the four fire sea¬ 
sons (Fig. 3: a is the whole fire season, b is the spring fire season, 
c is the summer fire season, d is the autumn fire season). It is 
found that good linear relationships exist between the two vari¬ 
ables for the four fire seasons with regression coefficients rang¬ 
ing from 0.16 to 0.61. It is interesting to note the linear relation¬ 
ship between area burned and the SSR value for the autumn fire 
season is the strongest among the four fire seasons, with the re¬ 
gression coefficient 0.61 as the highest. 
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Fig. 3 The relationship between area burned and mean SSR values 
during 1961-1990 period (a: the whole fire season; b: spring fire sea¬ 
son); c: summer fire season; d: autumn fire season). 

Time series of SSR over the Chinese boreal forest 

SSR was calculated for all eighteen meteorological stations, and 
these SSRs were interpolated to generate a SSR grid to cover the 
study region. Fig. 4 includes such interpolated grids for the mean 
30-year SSR values under both A2a and B2a scenarios, respec¬ 
tively. The SSR values ranged from 0.01-1.20 across the study 
region, and obviously SSR values showed an increasing pattern 
in the future compared to the baseline period. The fastest in¬ 
crease in the mean 30-year SSR happens in the summer fire 
season, and the second fastest increase is in the autumn fire sea¬ 
son. 


Conclusion 

In this study, we estimated the area burned in China’s boreal 
forest for the 21st century. Our results indicate that area burned 
would increase in the study region with continuing global 
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warming, and the results are similar to those found in previous 
studies conducted in Canada, USA, Russia, and Portugal. A sig¬ 
nificant relationship was found between historical area burned 
and the mean SSR values in the study region, which can be illus¬ 
trated by a simple linear regression with regression coefficient 
ranging from 0.16 to 0.61. Considerable evidences suggest that 
the increasing SSR value could translate into increasing area 
burned, and this finding is in accordance with that as found in 
Canada's boreal forest using two GCMs (Flannigan et al. 2005). 
According to the Hadcm3 model, the area burned would increase 
from 20% to 174 % across the most study region for the period 
2010-2099 under both scenarios of A2a and B2a. Furthermore, 
area burned would peak in 2080s under both scenarios, and area 


burned would increase by 20-174% under the A2a scenario and 
by 25-105% under the B2a scenario. The study region would 
experience 20-85% increase in area burned in the spring fire 
season under the A2a scenario, and 27-65% under the B2a sce¬ 
nario. However, the highest increase in area burned would ap¬ 
pear in the summer fire season. The study region may experience 
an increase in area burned of 174% under A2a scenario and 
106% under B2a scenario for 2080s. Alternatively, the region 
would face an increase of 22-174% in area burned in the spring 
fire season under the A2a scenario, and an increase of 25-106% 
under the B2a scenario. The results provided in this study may 
provide important information on how area burned is influenced 
by climatic changes for the fire management in the study region. 
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Fig. 4 Seasonal mean SSR (SMS) values for the baseline and the three future under A2a scenarios (a: the whole SMS for the baseline period; b: 
spring SMS for the baseline period; c: summer SMS for the baseline period; d: autumn SMS for the baseline period; e: the whole SMS for 2020s; f: 
spring SMS for 2020s; g: summer SMS for the 2020s; h: autumn SMS for the 2020s; i: the whole SMS for 2050s; j: spring SMS for 2050s; k: summer 
SMS for the 2050s; 1: autumn SMS for the 2050s; m: the whole SMS for 2080s; n: spring SMS for 2080s; o: summer SMS for the 2080s; p: autumn SMS 
for the 2080s). 
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